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Abstract 

We search for inclusive high Et diphoton events with large missing transverse 
energy in pp collisions at y / s=1.8 TeV. Such events are expected from pair 
production of charginos and neutralinos within the framework of the minimal 
supersymmetric standard model with a light gravitino. No excess of events is 
observed. In that model, and assuming gaugino mass unification at the GUT 
scale, we obtain a 95% CL exclusion region in the supersymmetry parameter 
space and lower mass bounds of 150 GeV/c 2 for the lightest chargino and 
75 GeV/c 2 for the lightest neutralino. 
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Supersymmetric models with a light gravitino (G), first proposed by Fayet |l|], have gen- 
erated recent theoretical interest . These models are characterized by a supersymmetry 
breaking scale A as low as 100 TeV and a gravitino which is naturally the lightest supersym- 
metric particle (LSP). The lightest superpartner of a standard model particle, assumed here 
and in most analyses to be the lightest neutralino (xT), is the next-to-lightest supersymmet- 
ric particle (NLSP). If Xi has a non-zero photino component, it is unstable and decays into 
a photon plus a gravitino (x? — > jG). 

In this Letter, we present a direct search for supersymmetry with a light gravitino in 
the framework of the minimal supersymmetric standard model (MSSM). In this framework 
the gaugino-Higgsino sector (excluding gluinos) is described by four parameters: Mi, Mi, 
fi, and tan/?, where Mi and M 2 are the U(l) and SU{2) gaugino mass parameters, /z is 
the Higgsino mass parameter, and tan/3 is the ratio of the vacuum expectation values of 
the two Higgs doublets ||. With the assumption of gaugino mass unification at the GUT 
scale, Mi = |M 2 tan 2 ^vK! where By? is the weak mixing angle. There are four neutralinos 
(x^, i = 1,2,3,4) and two charginos (xf, j = 1,2) whose masses and couplings are fixed 
by M%, fJ, and tan (3. We assume tan/3 > 1 in this analysis. 

We search for neutralino and chargino pair production in y/s=1.8 TeV pp collisions at 
the Fermilab Tevatron. The Xi is assumed to be short-lived, decaying within the detector to 
jG with a branching ratio of 100%. Decay to a Higgs boson is assumed to be kinematically 
inaccessible. R-parity conservation is assumed so that supersymmetric particles are pair 
produced and the LSP is stable and non-interacting. Thus pair production of charginos and 
neutralinos yields 77-^r events with high transverse energy (Et) photons and large missing 
transverse energy ($r), with or without jets. The high Et photons and large IpT provide a 
powerful tool for suppressing backgrounds. 

Recently D0 reported a search || for 77-^r events based on supersymmetry models with 
Xi as the LSP. In this analysis, we present the first experimental study of pp — > 77^t + X 
based on the MSSM with a light gravitino as the LSP. Using this model and more efficient 
photon identification and event selection criteria than in Ref. 0, we set the strongest limits 
to date in the supersymmetry parameter space, exceeding those from LEP experiments ||. 

The data used in this analysis were collected with the D0 detector during the 1992-1996 
Tevatron run at y/s=1.8 TeV and represent an integrated luminosity of 106.3 ± 5.6 pb _1 . A 
detailed description of the D0 detector can be found in Ref. [[7] . The trigger requires one elec- 
tromagnetic (EM) cluster with transverse energy Et > 15 GeV, one jet with Et > 10 GeV, 
and I$t> 14 GeV ($t> 10 GeV for about 10% of the data taken early in the Tevatron run). 
The jets in the trigger include non-leading EM clusters. Photons are identified through a 
two-step process: the selection of isolated EM energy clusters and the rejection of electrons. 
The EM clusters are selected from calorimeter energy clusters by requiring (i) at least 95% 
of the energy to be deposited in the EM section of the calorimeter, (ii) the transverse and 
longitudinal shower profiles to be consistent with those expected for an EM shower, and 
(iii) the energy in an annular isolation cone from radius 0.2 to 0.4 around the cluster in 
rj — 4> space to be less than 10% of the cluster energy, where r\ and 4> are the pseudorapidity 
and azimuthal angle. Electrons are removed by rejecting EM clusters which have either 
a reconstructed track or a large number of tracking chamber hits in a road between the 
calorimeter cluster and the event vertex. I$t is determined from the energy deposition in 
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the calorimeter for \rj\ < 4.5. 

To be selected as jjfa candidates, events are first required to have two identified pho- 
tons, one with > 20 GeV and the other with Ej? > 12 GeV, each with pseudorapidity 
|t7 7 | < 1.2 or 1.5 < |t7 7 | < 2.0. We denote the 28 events passing these photon requirements 
as the 77 sample. We then require fa> 25 GeV with at least one reconstructed vertex in 
the event to ensure good measurement of fa. No requirement on jets is made. Two events 
satisfy all requirements. 

The principal backgrounds are multijet, direct photon, W + 7, W + jets, Z — ► ee, and 
Z — > tt — > ee events from Standard Model processes with misidentified photons and/or 
mismeasured fa. The background due to fa mismeasurement is estimated using events 
with two EM-like clusters which satisfy looser EM cluster requirements than those discussed 
above, and for which at least one of the two fails the EM shower profile consistency require- 
ment (ii) above. In addition, these events must pass the photon kinematic requirements. 
These events, called the QCD sample, are similar to those of the 77 sample and are expected 
to have similar fa resolution. By normalizing the number of events with fa < 20 GeV in 
the QCD sample to that in the 77 sample, we obtain a background of 2.1 ± 0.9 events due 
to fa mismeasurement for fa > 25 GeV. 

Other backgrounds are due to events with genuine fa such as those from W+'^ (where 
'7' can be a real or a fake photon), Z — > tt — ► ee, and tt — > ee + jets production. These 
events would fake j'jfa events if the electrons were misidentified as photons. We estimate 
their contribution using a sample of e'7' events passing the kinematic requirements, including 
that on fa. Electrons are selected from the identified EM clusters with matched tracks. 
Taking into account the probability (0.0045) that an electron is misidentified as a photon, we 
estimate a background of 0.2±0.1 events. Adding the two background contributions together 
yields 2.3 ± 0.9 events. The fa distributions of the 77 sample and the background sample 
are compared in Fig. [I]. Also shown are the expected distributions from supersymmetry for 
two representative points in the (//., M2) parameter space. 

Chargino and neutralino pair production and decay are modeled using the spythia pro- 
gram ||, a supersymmetric extension of the pythia 5.7 program ||. Squarks and sleptons 
are assumed to be heavy. This assumption is conservative because light sleptons would lead 
to events with less jet activity and would therefore improve detection efficiency For light 
squarks, no change in efficiency is expected. To explore the parameter space, we choose to 
work in the (/1, M 2 ) plane while keeping tan/3 fixed. We generate XiXji xlxf an d xtxf 
events for a large number of points in the (//, M 2 ) parameter space. Table | shows the re- 
sulting theoretical cross sections a t h f° r several representative points, calculated using the 
CTEQ3L parton distribution function |T0j . To determine the signal efficiencies, Monte Carlo 



events are run through a GEANT |TT| based D0 detector simulation program, a trigger simu- 
lator, and the same trigger requirements, reconstruction, and analysis as the data. The total 
signal efficiency e (including efficiencies of the trigger, reconstruction, photon identification, 
and kinematic requirements) varies greatly, from ~ 0.01% to ~ 26%, depending largely on 
the masses of xt and xl an d their mass difference. For large masses such as those in Table |, 
the total efficiency e > 15%. The estimated systematic error on the total efficiency is 0.06e. 

With two events observed and 2.3 ± 0.9 events expected from background, we observe 
no excess of events. We compute 95% CL upper limits on the cross section a for the Monte 
Carlo sampled points in the (fj,, M 2 ) plane using a Bayesian approach [|I2[] with a flat prior 
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FIG. 1. The I$t distributions of the 77 and background samples. The number of events with 
lfir< 20 GeV in the background sample is normalized to that in the 77 sample. Also shown are 
the expected distributions (multiplied by 10) from two representative points in the supersymmetry 
parameter space, with tan (3 = 2. 



distribution for the signal cross section. The calculation takes into account the errors on the 
luminosity, the efficiency, and the number of background events. Depending on the values 
of the supersymmetry parameters, the 95% CL upper limits on the total cross section vary 
widely from several hundred pb for light charginos/neutralinos to a ~ 0.18 pb for heavy 
charginos/neutralinos. The upper limit <7d quoted in Table | is for events satisfying the 
kinematic cuts of this analysis at the generator level; comparison of a and op indicates the 
fraction of events yielding detectable particles for the various parameter points. 

To derive bounds in the (/x, M 2 ) plane, the values of \i and M 2 are varied around the 
sampled points until the theoretical cross sections a t h exceed the upper limits a. The 
interpolated bounds in the (/i, M 2 ) plane are shown in Fig. [2] for tan (3 = 1.05,2, 100. The 
regions below the lines are excluded by this analysis. The bounds depend on the value of 
tan f3 slightly, becoming stronger in the \i < half-plane and weaker in the other half-plane 
as tan f3 is increased. 

Figure [| compares the bounds in the (fi, M 2 ) plane for tan/3 = 2 with those estimated 
from LEP data within the framework of a light gravitino and assuming a 75 GeV/c 2 
selectron for t-channel exchange. These bounds exclude the region of parameter space 
suggested in Ref. || for the chargino interpretation of an event candidate shown by the CDF 
Collaboration JIB] . Also shown are the contours of constant mass for m^± = 150 GeV/c 2 and 
m^o = 75 GeV/c 2 . Since these are the largest masses for which the mass contours lie entirely 
in the excluded region, we obtain 95% CL lower mass limits of 150 GeV/c 2 for the lightest 
chargino and 75 GeV/c 2 for the lightest neutralino. This 75 GeV/c 2 lower mass limit also 
rules out a large part of the parameter space suggested for the selectron interpretation of the 
CDF event candidate in the model, as discussed in Ref. || . These mass limits are insensitive 
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GeV 


M 2 
GeV 


GeV/c 2 


m_± 
*i 


(Pb) 


Efficiencies (%) 


Limits (pb) 
a a D 


-160 


300 


143.9 


167.8 


0.12 


26.0±1.4 


36.4 


0.18 0.13 


-600 


140 


72.5 


146.4 


0.36 


17.2±1.2 


32.1 


0.28 0.15 


-800 


165 


84.7 


170.0 


0.20 


15.1±1.1 


26.4 


0.32 0.18 


200 


300 


118.1 


160.2 


0.15 


21.3±1.3 


31.9 


0.23 0.15 


400 


190 


89.4 


166.4 


0.19 


20.1±1.3 


32.5 


0.24 0.15 


800 


170 


83.2 


161.6 


0.25 


19.6±1.3 


33.4 


0.25 0.14 



TABLE I. Representative points in the (fi, M 2 ) plane for tan/3 = 2 with GEANT simulation. 



These points are chosen to be near our 95% CL bounds, where the experimental 95% CL cross 
section a equals the theoretical cross section oth- The efficiency e is for observing the total cross 
section a, while cd and <jd are the efficiency and cross section for observing the detectable events, 
those which satisfy the kinematic cuts > 20 GeV, E]? > 12 GeV, |r/ 7 | < 1.2 or 1.5 < |r/ 7 | < 2.0, 
and tfr > 25 GeV at the generator level. The total efficiency e = eo x ex where ex is the efficiency 
of the kinematic cuts. 
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FIG. 2. 95% CL bounds in the (ji,M 2 ) plane for tan/? 
line), and tan/3 = 100 (dashed line). 
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FIG. 3. Bounds in the (/i, M%) plane for tan/3 = 2. The region below the two solid lines is 
excluded at 95% CL. Also shown are the bounds estimated in Ref. || from LEP data (dotted line) 
and the contours of constant m-± = 150 GeV/c 2 (dashed line) and m^o = 75 GeV/c 2 (dot-dashed 
line). The hatched areas are suggested in Ref. Q for the chargino interpretation of the CDF event 
candidate in the model. 



to the choice of tan j3, varying less than 2 GeV/c 2 over the range 1.05 < tan (3 < 100, as long 
as our assumption that is the NLSP is satisfied. For large tan (3 values, this assumption 
may not be satisfied 0]. 

Most of the theoretical cross section for the 77-^r process is due to X1X1 an d xfx® 
production. For the large part of the parameter space with |/i| > M2, the relation m-± m 
m^o rs 2 x m^o holds, so we can express our cross section limits simply in terms of m^±. 
Figure [| shows the 95% CL upper limits for both processes, together with the theoretical 
predictions for tan/3 = 2 and /i = —500 GeV. The experimental limits are insensitive to 
the choice of tan/3 and \i while the theoretical cross section varies by about 10%. Our data 
rule out chargino masses below m 137 GeV/c 2 in models with a light gravitino, assuming 

> M 2 . This limit, though weaker than the 150 GeV/c 2 limit determined above from 
all processes contributing to final states, is useful for comparison with semi-exclusive 

calculations of gaugino production. 

In summary, we have searched for inclusive high Et diphoton events with large missing 
transverse energy. Such events are expected in the framework of supersymmetric models 
with a light gravitino. No excess of events is found. The null result, interpreted in this 
framework and with the assumption of gaugino mass unification at the GUT scale, yields a 
95% CL lower mass limits of 150 GeV/c 2 for the lightest chargino and 75 GeV/c 2 for the 
lightest neutralino. 

We thank G.L. Kane for many useful discussions and S. Mrenna for his help in using the 
SPYTHIA program. We also thank the staff's at Fermilab and collaborating institutions for 
their contributions to this work, and acknowledge support from the Department of Energy 
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FIG. 4. Measured 95% CL upper limits and predicted theoretical cross sections for \\ Xi an d 
X^Xo production as a function of m~± , assuming m~± ~ m,-o ~ 2 x m A -o . The vertical hatched line is 
the 95% CL lower limit on m-± determined using the total cross section for all chargino/neutralino 
pair production and all possible \i values, as determined in this paper. 

and National Science Foundation (U.S.A.), Commissariat a L'Energie Atomique (France), 
State Committee for Science and Technology and Ministry for Atomic Energy (Russia), 
CNPq (Brazil), Departments of Atomic Energy and Science and Education (India), Col- 
ciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and 
CONICET and UBACyT (Argentina). 
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